Spatio-temporal anomalous diffusion in heterogeneous media by NMR 
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For the first time, the diffusion phase diagram in highly confined colloidal systems, predicted 
by Continuous Time Random Walk (CTRW), is experimentally obtained. Temporal and spatial 
fractional exponents, a and n, introduced within the framework of CTRW, are simultaneously 
measured by Pulse Field Gradient Nuclear Magnetic Resonance technique in samples of micro- 
beads dispersed in water. We find that a depends on the disorder degree of the system. Conversely, 
[s, depends on both bead sizes and magnetic susceptibility differences within samples. Our findings 
fully match the CTRW predictions. 
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Anomalous diffusion (AD), that predicts the mean- 
square displacement (MSD) of a diffusing particle to grow 
nonlinearly in time, ([x(t) — x(0)} 2 ) oc t v (with v ^ 1), 
is a property of many complex systems and its related 
phenomena have been observed in various physical fields 
[5J |H1 HH HH H7J . As shown by Metzler and coworkers 
[TU] , the features of AD can be described using fractional 
calculus. According to these Authors, molecular AD in 
media can be described by defining the motion propa- 
gator (MP) as the solution of fractional diffusion equa- 
tions, which arise from the continuous time random walk 
(CTRW) model. These equations involve two fractional 
exponents, a and /i, which are, the orders of the time 
and space fractional derivatives, respectively [15] . The 
theoretical framework of CTRW is well established and 
has been corroborated by huge amounts of Monte Carlo 
simulations (see for example [ZQ) together with several 
experimental evidences, mainly obtained by using fluo- 
rescent spectroscopy |13l 119] . However, to our knowl- 
edge, an experimental a vs [i phase diagram, showing 
the competition between superdiffusion and subdiffusion 
of diffusing molecules in heterogeneous media, has never 
been carried out. 

In this Letter, we experimentally challenge the Metzler ct 
al. prediction of the CTRW theory [10] by measuring the 
fractional exponents a and [i by means of Nuclear Mag- 
netic Resonance (NMR) methods, based on Pulse Field 
Gradient (PFG) technique [18] , and by providing their 
interplay for the first time. 

In the last few decades, NMR water diffusion measure- 
ments have been a topic of extensive research, having 
broad applications in biophysics [2 and medicine [TJ. 
Translational self-diffusion in liquid systems can be mea- 
sured using PFG techniques, by applying magnetic gra- 
dient pulses (named diffusion gradients) to the system, 
in addition to the static magnetic field of the instrument 
itself. The signal attenuation, that depends on both dif- 
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fusion gradient strength and diffusion time, is simply the 
Fourier Transform (FT) of the MP. When MP is Gaus- 
sian, NMR signal attenuation follows a mono-exponential 
Stejskal- Tanner decay [18]. Conversely, when the motion 
is described by a non-Gaussian propagator, the signal at- 
tenuation deviates from a mono-exponential decay [9"Ill2j. 
In the present study, we have measured the fractional 
exponent a collecting the PFG signal attenuation as 
a function of diffusion times compared to the asymp- 
totic expression of FT of the MP for subdiffusive regime, 
obtained from [10] , Moreover, [i was extracted from 
the PFG signal attenuation as a function of gradients 
strength fitted to the expression of FT of the MP for su- 
pcrdiffusive regime, obtained again from [10] . PFG NMR 
experiments were performed on controlled porous media 
characterized by packed polystyrene micro-beads (Mi- 
crobeads AS, Norway) of various sizes dispersed in wa- 
ter. Mono-dispersed (mono-d) and poly-dispersed (poly- 
d) samples were used to check the potential ability of 
a and [i parameters to detect different degrees of sys- 
tem disorder. Moreover, internal magnetic field gradient 
(Gint) measurements were performed to investigate the 
influence of the magnetic susceptibility differences (A% m ) 
on measured a and /i parameters. Experimental results 
demonstrate that the exponent a strongly depends on 
the disorder degree of the system, and is independent of 
A^m- Conversely, the exponent [i is strongly correlated 
with Axm and pore sizes. In summary, mono-d samples 
are characterized by both ordinary diffusion and spuri- 
ous magnetic susceptibility effect, while poly-d samples 
show a real subdiffusion of water molecules. Thanks to 
peculiar features of NMR technique, an experimental ev- 
idence of the spatio-temporal competition between long 
rests and long jumps of diffusing molecules, as defined in 
[TO] , is reported here. 

Theory Taking into account the CTRW model [TO] , 
each displacement in a diffusion process can be charac- 
terized by two parameters: a waiting time r elapsing be- 
tween two consecutive steps, and a variable displacement 
length £. r and £ are independent random variables, dis- 
tributed according to the probability densities ip{r) and 
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A(0- When V(r) ~ r^-* 7 , with < a < 1, and A(£) 
is characterized by a finite variance, the characteristic 
waiting time T = / °° Tip(r)dr diverges and the resulting 
CTRW is subdiffusive. As a consequence, the asymptotic 
behaviors of the FT of the MP is 
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where a = er and, if cr > 1, a = 1. Under these 
conditions, the MSD grows sublinearly in time, namely 
(x 2 (t)) ~ i a . Fice-versa, when A(£) ^ |^r 1 " p , with 
< p < 2, and "0( T ) is characterized by a finite mean 
value, the jump length variance £ 2 = J_ £ 2 A(£)ci£ di- 
verges and a superdiffusive (Levy flight) CTRW is ob- 
tained. In this case, the behavior of the FT of the MP 
is 



W(k,t) ~ e 
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where p = p and, if /? > 2, /i = 2. In the case a < 1 and 
/i < 2, when spatio-temporal coupling of jump length 
and waiting time is taken into consideration, an a vs p 
phase diagram can be drawn, as reported by jlOj . In the 
specific case shown here, i.e., mono-d and poly-d micro 
bead samples, we can expect to detect subdiffusion (i.e. 
characterized by a < 1 and p — 2) and Gaussian (i.e 
characterized by a = 1 and p = 2) diffusion processes 
only. 

In NMR spectroscopy, the effect of diffusion on an en- 
semble of N magnetized spins can be described in terms 
of the behavior of spin phases, [$(rj)]j=i...jv, subjected 
to an applied magnetic field gradient 0. The PFG NMR 
technique is characterized by radio frequency pulses to- 
gether with a couple of magnetic field gradient pulses of 
duration 5, strength |g|, and separated by a delay time 
t = A. The first gradient pulse effectively phase encodes 
the molecular spins according to each molecule position. 
The second decoding gradient pulse completely reverses 
the spin phase if and only if the molecule has not diffused 
during the time A. Conversely, if the molecule diffuses, 
the intensity of the signal due to the molecule will be 
attenuated because of a non complete phase reverse. In 
case of ordinary diffusion, the degree of attenuation is a 
function of both g and A, and occurs at a rate propor- 
tional to the self-diffusion coefficient D. One of the most 
interesting features of PFG is that the measured signal, 
i?A(k), is the FT of the averaged propagator, P S (R, A), 
i.e. i?A(k) oc W(\t, A). Since A defines the time win- 
dow in which the diffusion processes are observed, PFG 
NMR techniques allow the measure of effective a and 
p values (a e ff and p e ff) using the aforementioned ex- 
pressions ([!]) and ([2]), respectively. In particular, a e // 
is obtained by a fitting procedure of experimental data 
measured by PFG at different A values to the expression 
0. The coefficient Peff is instead obtained by a fitting 
procedure of experimental data at different k = ^JgS 



values, where 7 is the nuclear gyromagnetic ratio. Tak- 
ing into account the A values usually selected to perform 
PFG experiments, we expect to measure p e ff = 2 and 
a eff < 1) depending on the investigated heterogeneous 
sample. Another peculiar property of the NMR signal in 
heterogeneous systems is its dependence on Ax m gener- 
ated at the interface between regions with different Xm- 
The presence of different susceptibilities can be quanti- 
fied through G int measurements, which can be extracted 
by the signal decay, generated by a Spin Echo (SE) NMR 
sequence [HI E] ■ As a consequence, both a e f / and p e f / 
could be, in principle, affected by A% m . 

Methods and Materials All measurements were per- 
formed on a Bruker 9.4 T Avance system, operating with 
a micro-imaging probe (10 mm internal diameter bore) 
and equipped with a gradient unit characterized by a 
maximum gradient strength of 1200 mT/m, and a rise 
time of 100 ps. As a first step, relaxation-time NMR 
measurements were performed to fully characterize the 
investigated porous samples. In order to determine the 
porosity p, as previously reported pQ , a SE imaging ver- 
sion using MSME (Multi Slice Multi Echo) sequence (rep- 
etition time TR = 1500 ms, matrix 128 x 128, slice thick- 
ness STH = 1 mm, in plane pixel dimension 60 x 60 pm 2 , 
number of scan NS = 8) at various echo times TE (from 
2.8 ms to 300 ms) was used to obtain SE decay in dif- 
ferent regions of each sample. Conventional mean dif- 
fusion coefficient D of water in each sample was mea- 
sured by means of a spectroscopic Pulsed Field Gradient 
STimulated Echo (PGSTE, TE/TR = 18/3000 ms, dif- 
fusion gradient pulses delay A = 80 ms, diffusion gra- 
dient pulses duration 8 = 4.4 ms and diffusion gradient 
strength g applied along the x, y and z axes using 32 
gradient amplitude steps from 6 mT/m to 100 mT/m). 
Gi n t along x, y and z axes was measured by using a spec- 
troscopic SE sequence (TR = 1500 ms, NS — 8) with 
N = 64 data points (corresponding to 64 echoes refo- 
cusing every 2 ms from 1 ms to 125 ms), as previously 
described [4]. Mean value of internal magnetic field gra- 
dient (MGi n t) was obtained as the average value along 
x, y and z axes. The characteristic diffusion length £]j = 
(2D A) 1 / 2 and the characteristic length t = (^^;) 1/3 , 
as defined in |llj . were obtained from these preliminary 
measurements. In particular, £* establishes the spatial 
region around each bead, within which internal gradi- 
ents act, dephasing the spins. A spectroscopic PGSTE 
with 5 = 4.4 ms, g = 0.10 T/m (i.e. k = 22481 to" 1 ) 
along x, y and z axes, with TR — 2.5 s, NS = 32 and 48 
values of A in the range (0.020 -j- 1.0) s was used to col- 
lect data to be fitted to equation 0, in order to extract 
the ct e ff value along x, y and z axes. Vice-versa, a spec- 
troscopic PGSTE with A/5 = 80/4.4 ms, TR = 2.5 s, 
NS = 16 and 48 gradient amplitude steps from 0.026 
to 1.02 T/m along x, y and z axes was used to collect 
data to be fitted to the function |2]), thus obtaining a 
measure of the Me// value along x, y and z axes. The 
mean values of a e // (^ a e//) an d p e ff (Mp e ff) were 
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obtained by averaging: Ma eff = lJ2i=x,y,z( a eff)i and 

All fitting procedures were performed by means of 
Levenberg-Marquardt algorithm. 

Nine samples in total were carried out using polystyrene 
micro-beads with nominal average diameters of 30, 20, 
15, 10, 6 and 0.050 fim and characterized by a |A% m | = 
\Xm° ~ x^ lystyrene \ ~ 1-59 in SI units. Six 10 mm 
NMR tubes were filled up to a volume of approximately 
2 cm 3 with beads mono-d in a solution of polyoxyethyle- 
sorbitan-mono-laurat (Tween 20) 10 -6 M and deionized 
water (conductivity^ 10~ 6 mho /cm), and investigated 
1 month after their preparation. The sphere packing 
in each mono-d sample, was characterized by a mean 
"sphere density" -q — 0.72 ± 0.03. As a consequence, 
water diffused in porous systems comprised by intercon- 
nected pores of mean diameter £ s greater than the one 
of cannon-ball packing configuration, i.e. £ s > v/3 4 +1 d, 
where d is the bead mean diameter. Three other poly-d 
specimens were prepared in a similar way, using equal 
volume fractions of beads of varying sizes. The sam- 
ple called 140 + 80 + 40 + 10 + 6 fim, was made with 
micro-beads of 140, 80, 40, 10 and 6 fim, while two other 
samples called 140 + 40 + 6 fim (w) and 140 + 40 + 6 fim 
were made with micro-beads of 140, 40 and 6 fim inves- 
tigated 1 week (w) and 1 month (m) after their prepa- 
ration, respectively. The latter sample (rj = 0.41 ± 0.03, 
p = 0.59 ± 0.03) is the most disordered one, while the 
former {jf = 0.69 ± 0.03, p = 0.31 ± 0.03) is less disor- 
dered when compared to the 140 + 40 + 6 fim specimens, 
but more disordered when compared to mono-d samples. 
Finally, one tube filled up with free water was also used 
as control. 

Results and discussion The experimental version of 
the theoretical a versus fi graph suggested by Mctzlcr 
et al. 10J, for all ten samples investigated, is shown in 
FIG{T] Data points obtained from ordered mono-d micro- 
bead samples (black filled symbols) lie close to the dashed 
line a e ff = , indicating the border line between su- 
perdiffusion and subdiffusion regime zones. These points 
are identified by a similar Ma e ff value very close to 1, 
and by different values of Mfi e ff, ranging from 1.75 to 
2, depending on bead sizes. Specifically, data points as- 
sociated to 0.050 fim, 6 fim and 10 fim mono-d samples 
are localized in a subdiffusion regime zone, while 15 fim, 
20 fim and 30 fim ones lie in the superdiffusion regime 
zone, close to the border line. The distribution of these 
data matches the experimental results showed in FIG (2] 
where the dynamic feature of water in all mono-d sam- 
ples is described by the interplay between the character- 
istic lengths l D , £ s and I* [11. When £ s < £ D ,£* (i.e. 
bead size lower than 15 fim in the investigated samples), 
motional averaging regime occurs and spins explore the 
entire pore many times before a dephasing arises, so any 
local magnetic field variation is averaged out by diffusion. 
Conversely, when £* < £d,£s (i-e- bead size is equal to 
and higher than 15 fim in the investigated samples), the 
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FIG. 1: Ma e ff vs Mfj, e ff diagram for ordered samples 
(black filled symbols), disordered samples (empty symbols) 
and free water (star symbol). Dashed line 2Ma e ff — Mfi R ff 
represents the effective ordinary diffusion region in which 
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(x(t) 2 } oc t ^ — t. The two dashed lines Ma e ff = 1 and 
Mfi e ff = 2 delimit the two regions, representing the superdif- 
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FIG. 2: Characteristic lengths £d, is and £* as a function of 
micro-beads size. 



slow diffusion case occurs and spins experience a no to- 
tally averaged Gi nt , which produces an additional spin 
dephasing [111 . Vice-versa, data points obtained from 
disordered poly-d micro-bead samples (empty symbols) 
lie, as expected, in the subdiffusive region. Specifically, 
Ma e ff decreases as the degree of the disorder increases, 
while Mfi e ff does not discriminate between ordered and 
disordered samples, and assumes values slightly smaller 
than 2. The broader the distribution of spatial length 
scales which characterize poly-d samples, the broader the 
distribution of characteristic resting time scales. As a 
consequence, water diffusing in investigated poly-d sam- 
ples has a non negligible probability to experience long 
resting times and, thus, to show an anomalous subdiffu- 
sive behavior quantified by the decreasing Ma e ff value 
from the unity, as bead disorder increases. To un- 
derstand the physical meaning of the fi e ff parameter, 
which in principle accounts for superdiffusion processes, 
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FIG. 3: Mfj, ef f as a function of MG lnt , measured by PGSTE 
NMR sequence. Solid line is the regression line M/i^ff = 
(-0.1331 ± 0.0061) In(MGint) + (1.613 ± 0.016), R = 0.996. 
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FIG. 4: Ma e ff as a function of MG int , measured by PGSTE 
NMR sequence. Solid line is the regression line Ma e ff = 
(0.000 ± 0.010) MG in t + (0.967 ± 0.022), R = 0.051. 



that in the investigated micro-bead samples, the mea- 
sured fi e ff values less than 2 are only due to artifact 
effects generated by Ax m at the interface between beads 
and diffusing water, fndeed it is known that water be- 
havior in these kind of interconnected micro-pore systems 
is Brownian or subdiffusive [51150]. However, these arti- 
facts, which simulate superdiffusion, provide information 
about different pore sizes. These observations give a ba- 
sis for further investigations on the correlation between 
Me// and G int . 

Conclusion Here we have measured the characteris- 
tic parameters of diffusion phase diagram, a and fj,, by 
means of NMR in systems marked by the presence of 
spatio-temporal competition between long rests and long 
jumps of diffusing molecules. Following the basic com- 
mon knowledge on complex systems, we made ordered 
and disordered media using mono- and poly-dispersed 
micro-beads. The NMR measure of their disorder degree 
fully matches the predictions. Moreover, we highlighted 
that the real physical mechanism, which gives rise to an 
apparent superdiffusion, is due to Ax m between the solid 
and liquid phase in heterogeneous samples. Inferring 
disorder features in heterogeneous media from diffusion 
phase diagram obtained by NMR may be of paramount 
value in a variety of applications from oil-well logging and 
dynamics of polymers to the diagnosis and monitoring in 
vivo of many diseases in the human body. 



we guessed that the dependence of \x e f / on bead size may 
be caused by Axm at the interface between beads and dif- 
fusing water. To validate this assumption, we obtained 
the plot of FIG(3j which shows a strong linear correla- 
tion (R — 0.996) between M/j, e ff and ln(MGi nt ) for all 
mono-d samples. In particular, the higher the MGi n t, 
the lower the M/i e // value. Conversely Ma e ff does not 
depend on A% m quantified by MGi n t, as shown in FIG|4] 
In the light of these correlations, it is reasonable to as- 
sume that values of Mfj, e f / less than 2 are spuriously due 
to Axm, that introduces a pseudo-absorb/desorb process 
of water molecules at the interface between bead surface 
and water, instead of a real superdiffusion mechanism. 
This is also shown in FlG(2j which clearly shows that 
slow diffusion processes become predominant as pore di- 
ameter l s increases. In particular, local gradients impart 
a phase shift to the spins within £* region that adds up 
to the phase shift given by the diffusion gradient pulse. 

As a consequence, some spins contribute to increase 
the degree of PFG signal attenuation; differently, other 
spins, which can be in a very distant zone from the first 
ones, will acquire a phase which will help to increase 
the signal. Due to indistinguishable spins associated to 
water molecules, this scenario simulates a superdiffusion 
regime of water molecules whose signal disappears in one 
spot and appears in another. In other words, we suggest 
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